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We introduce a simple model of static networks, where nodes are located on a ring structure, 
and two accompanying dynamic rules of repeated averaging on periodic node states. We assume 
nodes can interact with neighbors, and will add long-range links randomly. The number of long- 
range links, E, controls structures of these networks, and we show that there exist many types of 
fixed points, when E is varied. When E is low, fixed points are mostly diverse states, in which 
node states are diversely populated; on the other hand, when E is high, fixed points tend to be 
dominated by converged states, in which node states converge to one value. Numerically, we observe 
properties of fixed points for various _B's, and also estimate points of the transition from diverse 
states to converged states for four different cases. This kind of simple network models will help us 
understand how diversities that we encounter in many systems of complex networks are sustained, 
even when mechanisms of averaging are at work, and when they break down if more long-range 
connections are added. 

PACS numbers: 89.75.Fb,89.75.Hc,87.23.Ge,05.45.Xt 



I. INTRODUCTION 

Complex networks have been a popular subject of sta- 
tistical physics for more than a decade, because of their 
versatility and rich dynamic behavior[Tl [2]. In addition, 
they might be the only way to explain some of the univer- 
sal features of various systems that can be represented by 
networks. Even though structures of static networks can 
be a research subject of its own, studies on dynamic net- 
works are also flourishing. Three approaches can used on 
models of network dynamics, depending on the charac- 
teristics of systems in question: (i) node dynamics, focus- 
ing on dynamics of node states, while link structures are 
static, (ii) link dynamics, focusing on dynamics of links 
and their structures, while nodes are stateless, or have 
static states, and (iii) coevolving dynamics, where both 
links and nodes evolve, possibly influencing each other. 
Note that when some entities are assumed to be static, 
it means that the characteristic times for their changes 
are considerably longer than those of dynamic entities. 
One field that has attracted many physicists recently as 
an application of dynamic complex networks is a study of 
social systems[3, ¥ , where nodes are humans or groups of 
humans, and their relationships define links, forming so- 
cial networks. Social entities infiuence each other, mostly 
in attracting ways, but one of the fundamental charac- 
teristics of our societies is the sustained diversities. How 
are they being sustained when social entities have ten- 
dencies to assimilate with each other? How and when 
do these diversities break down? These are the questions 
we try to answer using complex networks on problems of 
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opinion dynamics [SHU], cultural dynamics [HHS] and so 
on. 

Here we introduce a model of static networks and in- 
vestigate node dynamics on these networks using nu- 
merical simulations mostly, while we have already stud- 
ied the same dynamic rule on dynamic networks [TTl [T5] . 
This model has nodes with periodic states with repeated- 
averaging rules [51 [T^. and there are two types of links: 
short-range and long-range links. Long-range interac- 
tions are rapidly increasing in our societies, as we are 
recently seeing with ever-advancing technologies, and we 
will use the number of long- range links, E, as the control 
parameter of the model. In our model, E changes the net- 
work structure, and has an important role in determin- 
ing what types of fixed points are reached from random 
initial conditions. When there are not many long-ragne 
links, diverse states prevail, but diversities will disap- 
pear eventually and more converged states will appear. 
As we vary E, we will find several types of fixed points, 
and estimate the likelihood of spontaneous emergence for 
given types of fixed points, using Monte-Carlo simula- 
tions. Network models with periodic nodes can be ap- 
plied not only to social systems, but also to other domains 
as in synchronization problems in oscillator networks [17]. 
For example, there have been many works on dynamics 
of Kuramoto model[TSl [TS] on complex networks [2D1I26) . 
and one dynamic rule in our model is similar to that of 
Kuramoto model. 

The plan of this paper is as follows. In Sec. II, the net- 
work model and two dynamic rules are introduced, and 
dynamics of networks with a small number of nodes is 
examined. In Sec. Ill, numerical results are presented. 
We find different types of fixed points using Monte-Carlo 
simulations, and observe their properties as E varies. Fi- 
nally, in Sec. IV, we conclude with discussions. 
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FIG. 1. (a) TV nodes in an ID periodic lattice with only 
nearest-neighbor connections (d = 1 here), (b) Examples of 
A'^ = 4 and 5 when d = 1. We can add E long-range links, 
where i?jv,d is the maximum number of long-range connections 
for given A'^ and d. 
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FIG. 2. (Color online) For iV = 1000, we observe (a) the 
clustering coefficient, C(e), and (b) the average shortest path 
length, I/(e), versus e for 1 < d < 5 (averaged over 1000 
random network configurations). 



II. MODEL 



We describe a network model with simple examples in 
Fig. [ij In this model, the network has A'' nodes, and 
they are placed in a ring structure, a ID periodic lat- 
tice. Each node is connected to neighbors within dis- 
tance d, assuming N ^ d; then, there are dN links, 
called short-range links. Now we add E distinct links 
randomly, not between already-connected neighbors, and 
these links are long-range links (or shortcuts). The total 
number of links becomes dN -t- E, and the average con- 
nectivity k is 2{d + E/N). If we introduce the normalized 
parameter e = E/En^cL-, where E^^d is N{N — 1 — 2d)/2, 
the maximum number of distinct long-range links, this 
becomes a normalized control parameter of the network 
(0 < e < l)[27j. This model can be viewed as a modi- 
fied version of the small- world network by Newman and 
Watts [55]. When d = 0, networks will have the same 
structural properties as random networks while E is the 
total number of links. When e = and d > 1, the net- 
work is a regular network, and, as e increases, the net- 
work undergoes several structural phases including small- 
world networks for low e. 

In Fig. [21 we observe the clustering coefficient C(e) and 
the average shortest path length L(e) as we vary e, for 
the cases for TV = 1000 and 1 < d < 5. When e = and 



d > 1, we can easily find C(0) and L(0) as below, 
3(d-l) 



C(0) 



2(2d-l)' 



, , , q + l\ { N -5 \ N 



(1) 

(2) 



where q and r are the quotient and the remainder when 
dividing {N — 2 + 5)/2 by d, and 6 is defined as 



_ / (TV: even), 
1 (TV: odd). 



(3) 



For d > 1, C{e) does not change much until e reaches 2 x 
lO"'^ {E ~ 100), while L{e) is reduced to about L(0)/10. 
As e increases further, C(e) decreases sharply to almost 
and starts to increase at e ~ 2 x 10^^ {E ~ 10''), 
and finally reaches the fully-connected network at e = 1. 
On the other hand, L{e) monotonically decreases until it 
reaches 1 at e = 1. 

Now we introduce a state variable (j)i to each node i, 
and assume (j)i is periodic with the range of < (/)i < 1. 
Then the state space of all TV nodes, which can be repre- 
sented by an n-tuple, {(J)q, (jti, 02j ■ • • j (t>N-i), becomes an 
TV-dimensional torus (TV-torus). Now we impose a syn- 
chronous and deterministic dynamic rule with discrete 
time t {t = 0,1,2, . . .), and the rule can be represented 
by a map for each node i at each time step as below. 
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where cr is a coupling constant (0 < ct < 1, but we will 
assume a — 1 for simplicity here), Kij is an element of 
the adjacency matrix (1 when there exists a link between 
i and j, and otherwise), ki is the degree of node i, 
and A{(j)j — (f>i) represents the difference between (f>j and 
(f>i. Then, the next value is the average of all states of 
neighbors of i and i itself when cr = 1. {li a < 1, (jji 
will move closer to the average, and, if cr > 1, which is 
not allowed here, the next value of cj)i will overshoot the 
average.) 

If the state variable is not periodic, we can simply use 
A.{x) = X. But the state variable is periodic here, and 
the definition of averaging can become ambiguous. For 
example, if we have two values, 0.1 and 0.9, there are two 
middle points, and 0.5, and which is the average"! The 
answer is both can be right. If we calculate the average 
from the point of view at, let's say, 0.45, the average is 
0.5 because distances to two points from 0.45 are -0.35 
and 0.45, respectively. But from the point of view at 0.1, 
the average becomes 0, because distances to two points 
are and -0.2. Then, an average in view of one value 
can be different from those in view of other values. This 
means that we need an anchor value to find an average, 
and that averages are relative. In Eq. (|4]), each node i is 
using its state value (fn as the anchor, and is finding its 
own average at each time step. But there still exists an 
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node updates its state with the averaged value especially 
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FIG. 3. The function A{(j}j — (f>i) representing tiie difference 
between 4>j and (f>i. (a) Rule A with a threshold in Eq. ([5|. 
(b) Rule B in Eq. Q. The dashed line is sin[27r((^j-0i)]/(2^. 
as in the Kuramoto model. 



ambiguity when distances between two values are exactly 
0.5; for example, if the anchor value is and a value is at 
0.5, the distance can be either 0.5 or -0.5. But, if we set 
A(a;) = at x = ±0.5, this ambiguity can be avoided. In 
addition, A(a;) should be periodic [A(a; ± 1) = A(a;) for 
— 1 < x < 1], because a; is a periodic variable. 

One way to define A(a;) is as below [see Fig. [SJJa)], 



^a{x) 




(if x< -1 + 61), 
(if - 1 + 61 < a; < 
(if -6 <x<e), 
(if 6'<x< 1-61), 
(if 1-61 < a;), 



(5) 



where is the threshold for interaction (0 < < 0.5), 
and we will call this rule as Rule A. Thresholds have been 
used in models of hounded confidence^^^ [25] in opinion 
dynamics, and can influence the dynamics greatly. In this 
work, we set 9 as 0.5, the maximum value, for simplicity. 
Now we can find the next value of (j)i unambiguously. 
What happens if we repeatedly find averages using Eq. 
Q for a given initial condition for N nodes? Values will 
eventually converge to a fixed point in the state space, 
where the average among neighbors for each node is the 
same as the state value of itself. This rule has been used 
in previous works by the author [TTJ [15]. Note that, due 
to the discontinuities in Ayi(x), a slight change of the 
state for one node can induce abrupt changes of other 
nodes at the next time step. 

We can avoid these discontinuities, and introduce an- 
other form of A(a;) by modifying Eq. (15]) as shown in the 
next equation [see Fig. [SJJ' 



^b{x) 



x + l (if x< -0.75), 

-X - 0.5 (if - 0.75 <x< -0.25), 

X (if - 0.25 < a; < 0.25), 

-a; + 0.5 (if 0.25 <a;< 0.75), 

x-l (if 0.75 <x). 



(6) 



and we will call this rule as Rule B. This rule resembles 
the interaction term in the Kuramoto model [TSl, and we 
can use another form of A(x) as below. 



^k{x) 



1 

2^ 



sin(27rx). 



(7) 



where results will be similar to those from Rule B (will 
not be shown here). Under both rules, A and B, each 



when \(j)j 



< 0.25, and we can call both of them as 



the rules of repeated averaging of periodic variables. One 
difference is that |Ab| decreases as |x| increases from 0.25 
to 0.5 unlike Rule A, making the meaning of averaging a 
little different. This can be interpreted as another repre- 
sentation of the concept of bounded confidence, because 
the influence between nodes starts to decrease when the 
difference is greater than 0.25. Later we will see that this 
fact can also play an important role for stability of some 
fixed points. 

We know that the converged states (or consensus if 
we interpret (jji as an opinion of node i, or synchronized 
states if we interpret (pi as a phase of one of identical 
oscillators), where all nodes have the same value, are al- 
ways stable fixed points under both rules. But there are 
other kinds of fixed points, and one type of fixed points 
can be represented by 



bo ± i— (mod 1) for < i < TV, 



(8) 



where n is a positive integer (when n = 0, it becomes the 
converged state). We can call these states as period-n 
states, and n is a period number. If '+' ('— ') is used 
in Eq. ([s]), (f>i increases (decreases) with the same rate 
as i increases. Here we look at the stabilities of these 
period-n states for small A^'s. 

1. When N ^ 3, E ^ and d = 1 (fully connected), 
period-1 states [e.g., (0, |, |)] are fixed points under 
both rules; but they are stable under Rule A, while 
unstable under Rule B, because, if the difference 
between any pair of infiuencing nodes is greater 
than or equal to 0.25, where the Ab'{x) changes 
the sign, the state cannot be stable under Rule B. 

2. When N = 4, E = and d = 1, period-1 states 
[e.g., (0, |, |, |)] are fixed points under both rules; 
but they are stable under Rule A, but unstable un- 
der Rule B. 

3. When A^ = 4, S = 2 and d = 1 (fully connected), 
period-1 states are unstable fixed points under both 
rules. They are unstable under Rule A, because, if 
the difference between any pair of influencing nodes 
is equal to 0.5, where A^(a;) changes the sign when 
9 — 0.5, the state cannot be stable. No period-n 
state (n > 0) can be stable under Rule B when the 
network is fully connected, because there always 
exist a pair of nodes with the difference greater than 
0.25. 

4. When N = 5, E = and d = 1, period-1 states are 
stable fixed points under both rules. 

5. When N — 5, E — 5 and d — 1 (fully connected), 
period-1 states are fixed points under both rules; 
but they are stable under Rule A, while unstable 
under Rule B. 



6. When N = 10, E = and d = 1, both period-1 and 
period-2 states are stable fixed points under both 
rules. 

7. When A^ = 10, -B = and d = 2, period-1 states 
are stable fixed point under both rules; but period- 
2 states are stable fixed points under Rule A, while 
become unstable under Rule B, because there exist 
pairs of nodes with the differences at 0.4. 

We can easily see that, for fully-connected networks, con- 
verged states are the only stable fixed points under Rule 
B, while other stable fixed points can exist under Rule A. 
On the other hand, when E = 0, there exist stable fixed 
points other than converged states under both rules, but 
some stable fixed points under Rule A are unstable un- 
der Rule B. From above examples, it is easy to find that, 
when n > Y (Rule A) or n > ^ (Rule B), period-n 
states when E = become unstable. In the next section, 
we will find properties of all possible stable fixed points in 
more detail using numerical simulations, when N = 1000 
and d <ti N with various e values. 



III. NUMERICAL RESULTS 

Once an initial state is given, the state will converge 
to a fixed point eventually. Here we will categorize fixed 
points, and observe what types of fixed points are reached 
for given e under both Rule A and Rule B. Then, using 
the Monte-Carlo simulations, we will also estimate ratios 
of initial states, for which a certain types of fixed points 
will be reached. In other words, we will estimate average 
volumes of the basins of attraction for some types of fixed 
points, which form attractors. For example, if a network 
structure is given, the set of all converged states can be 
viewed as an attractor for the given dynamical rule, and 
states in a certain portion of the state space will reach 
converged states as time increases. If we choose Nr ran- 
dom initial conditions from uniformly distributed states 
and A^'o of them reach converged states, N^/Nr is the es- 
timated volume of the basin of attraction for converged 
states (the volume of the whole state space, iV-torus, is 
1). If No/Nr ~ 0, the volume of the basin of attraction 
is extremely small, and the converged state is unlikely 
to be reached spontaneously in the system even though 
the converged states are stable fixed points, unless initial 
conditions are chosen deliberately. On the other hand, if 
Ng/Nr ~ 1, the basin of attraction for converged states 
can be close to the whole state space. 

In our simulations, if E is either or E'^v.d (e = 
or 1), there is only one possibe network configuration. 
But if E is between and EN,d (0 < e < 1), there are 
many possible configurations; hence we pick both a net- 
work configuration for given E and an initial condition, 
randomly out of uniformly distributed choices, for each 
run. Then, No/Nr will represent estimates of the aver- 
age volumes of those basins of attraction for converged 
states over all possible network configurations. The same 
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FIG. 4. (Color online) For iV = 1000 and E = 0, we observe 
(a) examples of period-n states; (b) examples of typical fixed 
points when d > 1 for Rule A; and (c) ratios of period-n 
states out of 1000 random initial conditions under both Rule 
A (d = 1) and Rule B (various d's). 



method can be used on period-n states, too. We will still 
use Nr = 1000 for all our simulations. 

First, we find fixed points when N = 1000 and E = 
(see Fig. [4]). Typical fixed points here are period-n states 
defined in Eq. (|8|. For d = I under Rule A and any 
d under Rule B, period-n states are only fixed points 
found [see Fig. [ija)], but for d > 1 under Rule A, other 
types of fixed points exist as shown in Fig. W[h) [mean- 
ing that Eq. (Is]) is not the only solution that satisfy 
yi,(t>i{t-\- 1) = (^i(i)]. Actually they dominate because 
no period-n state was reached from 1000 randomly given 
initial conditions when A; > 1 under Rule A. These states 
have small locally- converaed\60\ segments with scattered 
values In Fig. p[c), we observe percentages of period-n 
states, including converged states (n = 0), out of 1000 
random initial conditions. The higher n, the less likely 
they will be reached from random initial conditions in 
most cases. Under Rule B, as d increases, the highest 
n that can be found in fixed points decreases, and the 
volume of the basin of attraction for converged states 
increases; for example, when d — 50, about 50% of ini- 
tial conditions reach converged states. For d — 2 and 3 
under Rule A, percentages of period-n states including 
converged states are all zero in our simulations. 

Next, we observe examples of fixed points under both 
rules as we vary e from to 1 for the case of A^ = 1000 
and d = 3 in Fig. [5] If we look at fixed points under 
Rule A, all fixed points are similar to states in Fig.Jll^b) 
when E is less than about 1000 as shown in Figs. pTa)- 
[sjc). As E increases further, gathered states like the one 
shown in Fig.jSJJd) start to appear, while diversely spread 
states still exist as in Fig.lsje) [both are at _B = 3162 (e ~ 
6.5x 10^^)]. Converged states starts to appear when E ~ 
4000, and starts to dominate at i? ~ 20 000 with about 
97% of all fixed points. At this point, another types of 
fixed points called grouping states start to appear. In 
these states, nodes are divided into separated groups of 
nodes with clustered state values, and they can be called 
as m-group states, where m stands for the number of 
groups. For the case of = 0.5 (Rule A), even-numbered 
groups are not reached [3T], while, under Rule B, grouping 
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FIG. 5. For A'" = 1000 and d = 3, we show examples of 
stable fixed points for various E values numerically. Network 
configurations and initial conditions are randomly picked for 
given E's. (a)-(li) Under Rule A, (i)-(l) Under Rule B. 



states are unstable, and will not be reached. When m = 
1, states are just converged states, so 3-group states are 
the ones that appear mostly [see Figs. I5jf) and Islg)]. 
Widths of state values for groups tend to decrease as E 
increases further, reaching at i? = En^^i, while grouping 
states with to > 5 also starts to appear at _E ~ 120 000 
[see Fig. [gh)]. 

Under Rule B, fixed points look a little different. First 
of all, stable fixed points are mostly locally converged 
throughout the whole structure, showing piecewise- linear 
shapes in the scatter plots in Figs. [5li), [slj), and Islk) 
when £' = 10 and 100. In these states, some can be 
gathered as in Fig. ^j) and others can be spread as in 
Fig. IsFk), even at the same E value. Under Rule B, 
converged states starts to appear at i? ~ 500, and reach 
100% at £: ~ 1200; while gathered states as in Fig. [sjl) 
a.t E = 631 and converged states coexist in this region. 
As E further increases, all states will lead to converged 
states. 

If we categorize fixed points based on our numerical 
experiments so far, they are as below. 

1. Converged states: All nodes have the same state 
value in these states, and they are stable in all cases 
even though the volume of the basin of attraction 
for these states can be almost zero when e <C 1. All 
converged states satisfy 



'JN-l- 



(9) 



A set of all possible converged states can be viewed 
as an attractor, and the volume of the basin of at- 
traction will be estimated numerically. 



2. Periodic states: In these states, (pi^s increase (or 
decrease) linearly as i increases along the struc- 
ture with n cycles, and they are called as period-7i 
states, where a positive integer n is in the range 
of < 71 < f (Rule A) or < n < ^ (Rule B). 
We exclude n — 0, because pcriod-0 states are con- 
verged state. These states are stable fixed points 
when E = (when there is no long-range link) . An 
equation for the period-n states is given in Eq. (|8| . 
For each n, there are two attractors, one with the 
set of increasing period-n states, and the other with 
the set of decreasing period-n states. Here we will 
estimate the sum of volumes of basins of attraction 
for these two attractors for each n. 

3. Grouping states: Nodes can be divided into sepa- 
rated groups of nodes with clustered state values, 
and they can be called as TO,-group states, where 
TO stands for the number of groups. We exclude 
TO = 1, because 1-group states are converged states. 
Grouping states can be found only when e ~ 1 un- 
der Rule A. If e = 1, where the network is fully 
connected, each group will consist of nodes with 
the same value. Groups don't have to be equally 
distanced; rather it depends on sizes of groups as 
shown in Fig. p[g) . 

4. Other states: There can be fixed points that don't 
belong to above categories, and nodes look mostly 
scattered in the scatter plots. We can find some 
fixed points when N is small. For example, in Fig. 
Iljb), cases with A^ = 4 and 5 are shown, and we 
can find some fixed points under Rule A. For the 
case of iV = 4, i? = 1 and fc = 1, (0, |, |, g) is a 
fixed point; for the case oi N ~ 5 E — 1 and k ~ 1, 
(0, jj, jj, jj, Yj) is a fixed point; and for the case 
of A^ = 5, £: = 2 and fc = 1, (0, i, i, i, |) is a 
fixed point. For A ^ 1 and E ^ I, these fixed 
points can be highly complicated (see Fig. Isl. If 
we categorize these states further, there can be two 
different types. 

• Gathered states: As shown in Figs. [SJJd) and 
[Sjj), state variables are gathered near one 
value, even though some states can be far from 
this value. 

• Spread states: As shown in Figs. [5je) and 
Isjk), state variables are spread well between 
and 1, even though we can find some patterns 
in the scatter plots. 

We can elaborate more on two types of other states 
mentioned above. Under Rule A, when E is less than 
about 5000, the whole spectrum of states between these 
two types of states can be reached. But for the approx- 
imate range of 5 x 10^ < i? < 1.5 x 10"' (10~^ < e < 
3 X 10~^), they are clearly separated (see Fig. 6). Under 
Rule B, though, states are locally converged throughout 
the whole structure; therefore states are represented by 



an one-dimensional curves in scatter plots. For low E^s, 
these curves are piecewise linear, and points where linear 
segments meet are the locations where long-range con- 
nections branch off. If E is greater, the numbers of these 
branching points increase, too, and the curves become 
more complicated. As in Rule A, gathered and spread 
states can coexist for the cases with the same E. 

To see how nodes are spread for these fixed points, we 
introduce a quantity s. First, we need an anchor point 
because we learned in the previous section that an anchor 
point is necessary to find the average value of (jji 's. Here 
we set the anchor point (ps as the point where the most 
nodes are gathered state-wise; in other words, if we find 
the probability distribution function, P{(j)), the anchor 
point (f)s is the point where P{4>s) is the maximum for a 
given fixed point. Second, we can find the average of all 
4>iS, {4>)s, viewed from 4>s, using /S.a{x) with 6 — 0.5, 

1 ^"^ 
(0). = 0, + - ^ A^(0j - 4>s) (mod 1). (10) 

Finally, s is defined as 
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j=o 



This is the standard deviation of node states viewed from 
(^s|31]. For converged states, s should be 0, while, for 
uniformly distributed states like period-n states, s is 1. 
For other states, if it is close to 0, the state is gathered 
near one value; and if it is close to 1, the state is diversely 
spread. (Note that 1 is not the maximum possible value.) 
Using this index, we can see distributions of various fixed 
points in more detail. 

In Fig. [6] we observe values of s using Monte-Carlo 
simulations, which show us what types of fixed points 
are more likely to appear spontaneously. Under Rule A 
[see Fig. [6la)], all fixed points are spread states when E 
is less than about 1000 (e ~ 2 x 10^'^). But a transition 
occurs in the approximate range of lO'^ < _E < 2 x 10'' 
(2 X 10^'^ < e < 4 X 10~^), and gathered states start to 
appear. At i? ~ 4 x lO"' (e ~ 8 x lO^'^), converged states 
start to appear, too. All three kinds of fixed points co- 
exist until E ~2 X 10'' (e ~ 4 x 10~^), where converged 
states are dominant (about 97%) and spread states no 
longer exist. At this point, grouping states start to ap- 
pear, and 3-group states are the ones that appear mostly. 
For 2 X 10" < £; < 5 X lO"* (4 x 10"^ < e < IQ-^), 
converged states dominate and reach about 99%. Gath- 
ered and spread states have all disappeared, and grouping 
states start to increase as E further increases. Grouping 
states with m > 5 also starts to appear at E ~ 1.2 x 10^ 
(e ~ 2.4 X 10~^). Grouping states are in the top-right 
corner of Fig. l6la). 

Under Rule B, gathered and spread states also appear 
when E is small, even though states are mostly one- 
dimensional in scatter plots as seen in Fig. [Sj Converged 
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FIG. 6. (Color online) We observe s for 1000 runs for various 
E values {E > 10) when A^ = 1000 and d = 3. For given 
E, both positions of long-range links and an initial condition 
were chosen randomly for each run. Each pixel represents 
a run, and a column of pixels, representing runs for given 
E, are ordered, so that the run with the lowest s is in the 
bottom. Blue regions represent converged states {s — 0), and 
red regions represent quite diverse states (s > 1). (a) Values 
of s under Rule A. (b) Values of s under Rule B. There are 
12 dots in graphs, indexed by alphabets, and they correspond 
to the states depicted in Fig. [S] 



states start to appear at -E ~ 500 (e ~ 10 ^), and their 
distributions reach 100% at E ~ 1200 (e ~ 2.4 x lO'^). 
Note that more than 10% of initial conditions reach con- 
verged states when E — a,s shown in Fig. |4]jc). But 
those converged states disappear before E reaches 10. 
When E is greater than 1200, only converged states are 
observed, meaning that the volume of the basin of attrac- 
tion for converged states is the volume of the whole state 
space. There are no grouping states here as discussed 
earlier. One interesting point is that, for the region of 
e for small-world behavior (when e is less than about 
10~^), reached fixed points are all diverse states without 
any converged state. 

Now we turn our attention to the case of d = 0, where 
all links are randomly attached. It is easy to see that a 
necessary condition for a network to have fixed points 
other than converged states is for the network to be 
cyclic. For acyclic networks, only converged states are 
fixed points. In our model for d > 0, short-range links 
guarantee that these networks have at least a cycle, and 
we observed many types of fixed points so far. What hap- 
pens if there is no prearranged short-range links for the 
case of d = 0? First, there can be more than one compo- 
nent, a group of connected nodes, for low i5's; hence we 
only consider states of nodes for the largest component 
when calculating s and Pi(t>), the probability distribution 
function of (/)i's for a given fixed point, since we are inter- 
ested in looking at how connected nodes influence each 
other. Second, the chance of the largest component of the 
random network to be acyclic is not negligible when E is 
low, but, as E increases, that chance becomes extremely 
small. But, for higher £"s, we expect to see similar re- 
sults as those for the cases of d > 0, because influences 
of short-range links become negligible when E ^ dN. 

In Fig. [71 we look at some examples of flxed points for 
the case of iV = 1000 and d = 0, and observe behavior of 
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FIG. 8. (Color online) We observe values of E at which 
converged states start to prevail under both rules, (a) Ratios 
of converged states out of 1000 runs versus k, the average 
connectivity, when d — and 3 under both rules. Values of k 
at which ratios reach 0.5 will be called as kc- (b) kc versus A*' 
for four cases. Slopes are exponents, ^'s, in kc oc A''^. 



FIG. 7. (Color online) We observe s for 1000 runs for various 
E values {E > 500) when TV = 1000 and d = 0. (a)-(d) Dis- 
tributions of states P{<j)) under Rule A. (e)-(h) Distributions 
of states P{4>) under Rule B. (i) Values of s under Rule A. (j) 
Values of s under Rule B. There are 8 dots in graphs, indexed 
by alphabets, and they correspond to the states in this figure. 



s as in Figs. [5] and [6J Since scatter plots of (pi versus i are 
no longer meaningful, we use the probability distribution 
function P{(j)) to represent each fixed point. In Figs. 
^a)-l7rh), some fixed points under both rules are shown. 
We can find all types of fixed points except the periodic 
states. We can also observe 2-group states under Rule A 
here, unlike the cases of d > [for 2-group states, a group 
usually consists of only few nodes as shown in Fig. mb)]. 
In Figs. IrFi) andl7|j), we observe s values for 1000 runs 
for each E under both rules. When E is less than about 
600 (e < 1.2 X 10^"^), converged states are dominant as 
stable fixed points, but, as E increases further, diverse 
states appear and later disappear with the same kind of 
transition from diverse states to converged states as we 
have already observed for cases of d = 3 in Fig. [6j 

One important dynamic property on networks we have 
been using is this transition from diverse states to con- 
verged states as the number of links E increases. (Un- 
der Rule A, grouping states appear after this transition.) 
Then, we can ask when this transition occurs. One way to 
quantify the answer is finding the point at which the ratio 
of converged states out of all runs becomes exactly a half, 
and we will use the average connectivity, instead of E, to 
represent this value: kc [note that k = 2{d + E/N)]. In 
Fig. Isla), we show ratios of converged states out of 1000 
runs versus k for four cases we have used at TV = 1000: 
d = 0, 3 (Rule A), and d = 0, 3 (Rule B). Under Rule A, 
two curves match very well, and kc is about 27; while un- 
der Rule B, two curves don't match because the numbers 
of nodes and links for the largest component are much 



smaller than N and i? at fc ~ fc^ when d = 0. (For the 
case of d = under Rule A, on the other hand, the pos- 
sibility of having more than one component is extremely 
low at fc ~ fee.) Finally, in Fig. IsFb), we observe how kc 
changes as N increases. Overall dynamic behavior will 
not be much different, but how kc varies with A'^? Numer- 
ical results show that kc fit well with N^, and exponents 
f 's for four cases are also found. Exponents don't seem 
to depend on d under Rule A, and ^ ~ 0.23; while under 
Rule B, ^ ~ 0.10 {d = 0) and ^ ~ 0.03 {d = 3). Resuhs 
suggest that more long-range links are needed to get con- 
verged states under Rule A, and kc increases faster with 
N, compared to results under Rule B. It is also interest- 
ing that, under Rule B (d = 3), /cc doesn't increase much 
with A'' (for example, kc ~ 7.4 at A/ = 10^ and kc ~ 7.8 
at Af = 10^*). 



IV. DISCUSSIONS 

Using a simple network model, we have found fixed 
points for periodic state variables using two dynamic 
rules of repeated averaging. Rule A and Rule B. This 
model has a parameter e (0 < e < 1), representing the 
number of long-range links, and stable fixed points and 
their likelihood of appearing spontaneously from random 
initial conditions have changed as e varies. We catego- 
rized fixed points as converged states, periodic states, 
grouping states, and other states (gathered and spread). 
Grouping states are only reached under Rule A, while 
other types of fixed points have been found under both 
rules. As the number of long-range links increases, the 
transition from diverse states to converged states oc- 
curred in all cases, and we represented the point of tran- 
sition using kc, the average connectivity when the ra- 
tio of converged states becomes exactly 0.5. Numerically 
found fee's fit well with A^^, and exponents, ^'s, have been 
found. For Rule A, we need more long-range connections 
(higher kc) to get converged states, and the exponent ^ is 



greater, compared to results under Rule B. In addition, 
diverse states can reappear in the form of grouping states 
under Rule A, when the network is close to being fully 
connected. 

The same dynamic rule (Rule A) has been used in 
previous works [TTJ [15] for dynamic interaction networks 
with random-walking nodes in two-dimensional lattices. 
In these networks, short-range links constantly change 
as nodes move randomly in the lattice, thereby giving 
stochasticity to the model. Unlike the static network 
used here, converged states are the only stable fixed 
points, and periodic and grouping states are found to be 
metastable with finite lifetimes. For dynamic cases, pe- 
riodic states are metastable when < e ^ 1, while they 
can be fixed points when e = here. In addition, fixed 
points like gathered and spread states cannot be stable in 
dynamic networks, because these states are only stable 
specifically for given network configurations. 

The motivation for this work is to answer the questions 
of how the diversities exist in reality even though many 
mechanisms of assimilation are at work, and how these di- 
versities disappear if we increase long-range interactions. 
For example, if we look at human cultures, they are still 



quite diverse, but neighboring cultures differ only slightly 
because neighbors usually interact and assimilate. One 
can use interaction networks to describe these systems, 
and argue that the prevalence of short-range interactions 
are the main cause of these kinds of diversities. Our 
model showed that, if there are only short-range links, 
periodic states, which are locally converged and yet di- 
verse, can be stable fixed points. The point of transition 
from diversity to convergence we found in this model can 
be used as a guide when we roughly determine how many 
long-range links are necessary to make a group of nodes 
converge to one value. Even though the reality is more 
complex, we hope this simple model can give us insight 
when we try to understand phenomena of sustained di- 
versities and their breakdown. 
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